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Abstract 

In this paper I introduce the probability distribution of the local overlap in spin glasses. The 
properties of the local overlaps are studied in details. These quantities are related to the recently 
proposed local version of the fluctuation dissipation relations: using the general principle of stochastic 
stability these local fluctuation dissipation relations can be proved in a way that is very similar to the 
usual proof of the fluctuation dissipation relations for intensive quantities. The local overlap and its 
probability distribution play a crucial role in this proof. Similar arguments can be used to prove that 
all sites in an aging experiment stay at the same effective temperature at the same time. 

1 Introduction 

Up to now in disordered systems the overlap and its distribution were considered as global quantities, that 
were defined for the whole system Q. However in systems with quenched disorder it is possible to define 
(in a non-trivial way) a local overlap that has a point-dependent probability distribution j2j. This new 
object has remarkable properties that we explore in this paper. 

One of the most interesting results is related to the local generalization of fluctuation dissipation relations 
(FDR) in off-equilibrium dynamics. 

It is well known that the fluctuation dissipation relations in off-equilibrium dynamics are a crucial tool 
to explore the landscape of a disordered system H ||, ||, [(J. These fluctuation dissipation relations are 
different from the prediction of the fluctuation dissipation theorem at equilibrium. They can be expressed 
in a rather simple form that can be easily interpreted from the theoretical point of view. Moreover the 
main parameter entering in the fluctuation dissipation relations has simple interpretation from the point 
of view of equilibrium statistical mechanics |^, ^, [l^] . This fact has consequence that for a given system 
the form of the FDR is universal in off-equilibrium dynamics, i.e. it does not depends on the details of the 
dynamics and on the way in which the system is put in an off-equilibrium situation (as soon the system 
remains slightly out of equilibrium). 

In the mostly studied case one considers observables that are the average over the whole sample |n], [j], 
P^4 ||[ |ll| [ll| . In this case the static-dynamic relations connect the FDR to the static average of global 
quantities. Recent there have been a few investigations on FDR that involve only given local variables 
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using the probability distribution of the single spin overlap it is possible to give 
a theoretical foundation to these local FDR and to derive the appropriate static dynamic relations, that 
involve the probability distribution of the local overlap Q . 

The paper is organized as follows: in the second section we define our main new theoretical tool: the 
probability distribution of the local overlap; we compute its properties in a few simple cases, in the third 
section we recall the usual global FDR, while in the next section we recall the proposed local FDR 0, |. 
In the fifth section we derive from general principle the local FDR, we prove the appropriate static-dynamic 
relations and show that in an ageing regime, in spite of the existence of local heterogeneities all sites at 
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a given time must be characterized by the same effective temperatures. In the next section, before the 
conclusions, we discuss some methods to compute the probability distribution of the local overlap. Finally, 
it the appendix, we present some consideration on systems with finite volume. 



2 The local overlap distribution 

2.1 Some heuristic considerations 

This section contains the definition of the local overlap probability distribution Pi (q) . 
Usually the overlap of two equilibrium configurations (<r and r) is defined as 

q = ^ L , (1) 
N y ' 

where qi = a^. For a given sample we can define the overlap probability distribution Pj(q), where J label 
the sample. In the glassy phase the function Pj(q) depends on the sample also for very large samples: the 
global overlap probability distribution Pj(q) is not a self-averaging quantity. The physically interesting 
quantity is 

P{q) = PM, (2) 

where the bar denotes the average over the couplings J. 

In a different approach JTo| one considers the response of the system to the appropriate perturbation 
and in this way one can define for a given sample a function P r (q), that for large systems should be self- 
averaging, i.e. J-independent. This new order parameter distribution codes the thermodynamic responses 
of the system to random perturbations. According to the principle of stochastic stability (t], [|, ^, |To) it 
should coincide with P(q). 

It is evident that the definition of the local overlap must be rather different from that of the global 
overlap. Indeed in a naive approach the local overlap of two equilibrium configurations (i.e. cr^rj) is always 
equal to ±1; if use a naive definition the probability distribution of the local overlap should be the sum of 
two delta functions at ±1: in this way one gets a trivial result. 

Moreover our aim is to define a site dependent, sample dependent Pi(q) such that 

PM= LI {Pi(9i)dqi)s(^--q) (3) 

i=l,N ^ ' 

is a self-averaging quantity. In other words P g {q) cannot coincide with Pj(q) but it should be equal to 
P r (q) . In the following we will extend the approach of |1(J in order to arrive to a definition of Pi (q) that 
satisfies the aforementioned properties. 



2.2 The definition of the local overlap 

Let us start from a spin glass sample and let us consider M identical copies of our sample: we introduce 
TV x M af variables where a = 1,M (eventually we send M to infinity) and N is the (large) size of our 
sample (i = 1, N). The Hamiltonian in this Gibbs ensemble is just given by 

Hk(<t)= ]T H(a a ) + eH R [a] , (4) 

a=l,M 

where H(a a ) is the Hamiltonian for a fixed choice of the couplings and the Hr[(t] is a random Hamiltonian 
that couples the different copies of the system. A possible choice is 

H R [a]= K »t + \ (5) 

a=l,M;i=l,N 

where the variables Kf are identically distributed independent random Gaussian variables with zero average 
and variance 1. In this way, if the original system was d dimensional, the new system has d+1 dimensional, 
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where the planes are randomly coupled. We can consider other ways to couple the systems (e.g. Hr[<j] 



T, a , b =i,M;i=i,N Ki'^fVi)- An other possibility 



is 



H R [a}= K '::>>': ( 6 ) 

a,b=l,M;i,j=l,N 

where the variables K are identically distributed independent random Gaussian variables with zero average 
and variance (NM) -1 . As we shall see later the form of Hr is not important: its task it to weakly couple 
the different planes that correspond to different copies of our original system. The first choice eq. (|J) is the 
simplest to visualize and it is the fastest for computer simulations, the last choice eq. (|^) is the simplest 
one to analyze from the theoretical point of view. In the following we do not need to assume a particular 
choice. 

Our central hypothesis is that all intensive self average quantities are smooth function of e for small 
e. This hypothesis is a kind of generalization of stochastic stability. According to this hypothesis the 
dynamical local correlation functions and the response functions will go uniformly in time to the values 
they have at e = 0. 

We now consider in the case of non-zero e two equilibrium configurations a and r and let us define for 
given K the site dependent overlap 

V rr a T a 

qi[a - T) = M (7) 

We define the K dependent probability distribution P^{q) as the probability distribution of the previous 
overlap. If we average over K at fixed e we can define 



P?(q)=Pf(q), (8) 
where the bar denotes the average over K. We finally define 

P i (q) = ]imP?(q) , (9) 



where the limit e — > is done after the limits M — > oo and N — > oo (alternatively we keep eM and eN 
much larger than 1). 

In order to be consistent we the usual approach we should have that if define 

It = — , (10) 

the probability distribution Pt(q) of qt should be self-averaging (i.e. J independent and it should coincide 
with the function P(q) that is the average over J of Pj(q): 

P t {q)=P g {q)^Ej[Pj{q)] . (11) 

This crucial relation will be proved in the next section. 

In the nutshell the construction is rather simple. We consider M weakly coupled copies of the original 
lattice and in the limit M going to infinity we can define local thermodynamics averages. We will assume 
that we stay also in the infinite volume limit. A discussion of what happens for finite N but M = oo will 
be presented in the appendix. 

2.3 Some heuristic considerations 

Which is the rational of this baroque construction? The heuristic idea is simple. For finite N and given J 
we can approximatively decompose the usual Boltzmann expectation value (•} into different states labeled 
by a: 

(■) = E^->« ( 12 ) 

a 

If we define the state dependent magnetizations 

m° = (ct^ . (13) 
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we have that the overlap among two states is 

Q , 7 = T,i=i,N m f m l (u) 

N K ' 

and that the usual overlap distribution is given by 

Pj(q) = Y,w a w 1 6(q a "-q) . (15) 

a, 7 

The average over J gives the usual P(q). 

However we do not want to average over the J;we want to stick to a given system. In order to make a 
sensible selfaveraging definition we have to consider an ensemble of system so we will consider the original 
system plus a small random perturbation H R . However, if the expectation value of the random perturbation 
in the various states (E R = (H R )/ a ) is much larger than one, the new w's will be proportional to 

w a exp(-/?£f ) (16) 

and the states will be completely reshuffled. In a first approximation we can assume that the set of states 
remains the same after the perturbation, but the thermodynamically most relevant states do change. The 
way in which they change depends on the distribution of the a. However the principle of stochastic stability 
tell us the we get the same function P(J) if we do the average over the random perturbation at fixed J or 
over the ensemble of the J in absence of the random perturbation. 

At non zero e we can assume that the systems we state is same state A with a probability wa and the 
state A is characterized by the fact that the a th system stay in the state a a and the states a are the same 
as at e = 0. By changing the variables K we change the weights wa- We can thus write 

EZ/o=l,M 1a a ,-y a , 1f7 v 
W A W C Jr ( 17 ) 

A,C 

It is clear that we have to prove that this way to generate the weights A and C is such the distribution of 
q t satisfies eq. eq. ( fill ) • This will be shown in the next section. 

2.4 Two explicit cases 

Let us consider here the previous construction in two cases where we can perform the relevant computations 
in an explicit way. 

The fist case we consider is one step replica symmetry breaking. Here the system may stay in states 
labeled by an index a. Each state is characterized by a free energy f a and the probability of finding a state 
in the interval [/, / + df] is given by 

N{f) = exp(m/3(/ - /„)) . (18) 

The states of the Gibbsian ensemble will be characterized by the weights wa- These states will have the 
same distribution of probability (eq. (jig)). However they will be characterized by variables a a that are 
different for each state. Therefore we can write: 

E—Cc a 7a 
a=l,M i i /ms 

«i = N (19) 

for A ^ C and 

.. . a=l,M °i °i , oa n 

Qi ' = jj . ( 20 ) 

where of is an equilibrium configuration of the state a. 

In the limit M — > oo we can use the law of large numbers and we find with probability one that 

A,C a,~i _ 

Qi =Qi -Qi , 

A. A a, a 1 fnt\ 

Qi =Q l =Qi , (21) 
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where the bar denotes the average over all the states, we have used the notation: 



qr = mfrn] (22) 
and mf the magnetization at the site i in the state a: 

mf = (<i t ) . (23) 

We finally find that 

Pi(q) = mS(q - g?) + (1 - m)S(q - q}) . (24) 

In other words the construction we have used performs automatically the average over all the possible states. 
It is remarkable that in this case the value of m is constant all over the system, and the site variability is 
present only in the values of q® and q\ . 

The same computations can be done in the case where replica symmetry is broken at two steps. In this 
case the states can be clustered into families, each state is labeled by two indices (e.g. c, 7) the first index 
label the family and the second index labels the state into the given family. In this case we have two free 



energy distributions of the form given in eq. (18) , one for the states (characterized by a parameter m s ) 
and the second for the families (characterized by a parameter mf), where m s > mf. In the same way as 
before we find 

PM = m f 5(q - q°) + (m s - m f )6(q - q]) + (1 - m s )5{q - q 2 ) , (25) 

where 

qr A5 = 9°, q^ S = q\, q^^ql- (26) 

We notice that in this approach we have two quantities: the weights, that are global quantities, and the 
magnetizations that depend on the point. For simplicity let us restrict our analysis to the one step replica 
symmetry breaking. In this case for each point we can reconstruct the probability distribution of qi if we 
know the probability distribution dPi(m) of finding a magnetization mi at the site i in a generic state: 

?° = ( I dPiim)™) , q}= fdP t (m)m 2 (27) 



It may be interesting to note that on the Bethe lattice, in the case where one step replica symmetry 
breaking is exact, the probability distribution Pi(m) depends only on the local environment (the coupling 
of the nearby points) and it can be computed in the large N limit by solving local equations |2^, p3[ , 
confirming the fact, that will be proved in full generality in the next section, that also the probability 
distribution Pi(q) depends on the local environment. 



3 The Global Fluctuation Dissipation Relations 

In this section we will find a short summary of the main results concerning the global FDR. 

The usual equilibrium fluctuation theorem can be formulated as follows. If we consider a pair of conju- 
gated variables (e.g. the magnetic field and the magnetization) the response function and the spontaneous 
fluctuations of the magnetization are deeply related. Indeed if R eq {t) is the integrated response (i.e. the 
variation of the magnetization at time t when we add a a magnetic field from time on) and C eq (t) is the 
correlation among the magnetization at time zero and at time t, we have that R eq (t) = [3(C eq (0) — C eq (t)), 
where (3 — (kT)^ 1 and 3/2 k is the Boltzmann-Drude constant. If we we eliminate the time and we plot 
parametrically R eq as function of C eq we have that 

<28) 

The previous relation can be considered as the definition of the temperature and it is a consequence of the 
zcroth law of the thermodynamics. 

The generalized fluctuation dissipation relations (FDR) can be formulated as follows in an aging system. 
Let us suppose that the system is carried from high temperature to low temperature at time and it is in 
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an aging regime. We can define a response function R(t w , t) as the variation of the magnetization at time 
t when we add a a magnetic field from time t w on, in a similar way C (t w , t) is the correlation among the 
magnetization at time t w and at time t. We can define a function Rt w (C) if we plot R(t w ,t) versus C(t w ,t) 
by eliminating the time t (in the region t > t w where the response function is different from zero. The FDR 
state that for large t w the function R tw (C) converge to a limiting function R(C). We can define 

-§=mC) (29) 

where X(C) — 1 for C > Coo = lim^oo C eq (t), and X(C) < 1 for C < Coo- The shape of the function 
X(C) give important information on the free energy landscape of the problem, as discussed at lengthy in 
the literature. 

It has been shown that in stochastically stable system the function X{C) is related to basic equilibrium 
properties of the system. Let us illustrate this point by considering for definitiveness the case of a spin 
glass. 

Spin glasses are characterized by the presence of a quenched disorder (i.e. the coupling J among spins). 
For a given probability distribution of the J, there are quantities that do not depend on the particular 
(generic) realization of the J and are called self-averaging: the response function and the correlation of the 
total magnetization belong to this category. On the contrary there are other quantities that depend on the 
choice of the J. A typical example of non self-averaging quantity is given by the probability distribution 
of the overlap. For fixed value of J, given two equilibrium configurations a and r, we define the global 
overlap as 

Q(v,t) = ^ (30) 

where TV is the total number of spins. Let us suppose that there is a magnetic field (albeit infinitesimal) 
in such a way that the overlap is positive, otherwise the overlap is usually defined as the absolute value of 
the previous expression. 

The probability of distribution of q is Pj(q) and it depends on J. The function P(q) is defined as the 
average of Pj(q) over the difference choices of the coupling J and obviously depends of the probability 
distribution of the variable J. 

It is convenient to introduce the function x(q) defined as 

x(q) = P P(q')dq' (31) 
Jo 

or equivalently 

Obviously x(q) = 1 in the region where q > qEA, where qsA is the maximum value of q where the P(q) is 
different from zero. 

The announced relation among the dynamic FDR and the statics quantities is simple 

X(C) = x(C) (33) 

We shall see later that this basic relation can be derived using the principle of stochastic stability that assert 
that the thermodynamic properties of the system do not change too much if we add a random perturbation 
to the Hamiltonian. All that is well known. In the nest section we shall se how to play the same music will 
local variables. 



4 The Local Fluctuation Dissipation Relations 

There are recent results that indicate that the FDR relation and the static-dynamics connection can be 
generalized to local variables in systems where a quenched disorder is present and aging is heterogeneous, 
We shall see that there finding need a more general framework to be explained. 
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For one given sample we can consider the local integrated response function Ri(t w ,t), that is the 
variation of the magnetization at time t when we add a magnetic field at the point i starting at the time 
t w . In a similar way the local correlation function d (t w , t) is defined to the correlation of among the spin 
at the point i at different times (t w and t). Quite often in system with quenched disorder aging is very 
heterogenous: the function d and Ri change dramatically from on point to an other. 

It has been observed in simulations that local fluctuation dissipation relations (LFDT) seems to hold 

-^=mo, (34) 

where -Xj(C) has quite strong variations with the site. 

It has also been suggested that in spite of these strong heterogeneity, if we define the effective /3?^ at 
time t at the site i as , 



dRi(t w ,t) 



= ^X i {t w ,t)=0f f {t w ,t) , (35) 



dCi (tyj , t) 

the quantity j3%^ does not depend on the site. In other words a thermometer coupled to a given site 
would measure (at a given time) the same temperature independently on the site: different sites are 
thcrmometrically indistinguishable. 

These empirical results calls for a theoretical explanation. The aim of this note is to show that these 
results are consequence of stochastic stability in an appropriate contest and that there is a local relation 
among static and dynamics. In the next section we will define in an appropriate way the local probability 
distribution of the overlap for a given system at point i (i.e. Pi{q)). We will also define the function Xi(q) 
as 

Xi(q) = f Pi(q')dq' (36) 
Jo 

and we will show that the static-dynamic connection for local variables is very similar to the one for global 
variables and it is given by: 

X i (C)=x i (C) (37) 

In order to arrive to to prove of the local FDR, that will be presented in the next section, is it convenient 
to reconsider the dynamic definition of the correlation function C, (t w ,t). It is clear that C{ (t w , t) cannot 
be measured by observing only one single history of our sample: <Ji{t w )(Ji{t w + 1) = ±1. The correlation 
function is obtained by averaging over all the possible history, i.e. by repeating the experiment M times 
and sending eventually M to infinity. In other words the two time local correlation function is not a self 
averaging quantity as far histories are concerned. 

If we want to define the correlation in such a way that it can be measured by observing a single history 
we have to consider the Gibbsian introduced in section II. We can consider M identical copies (or clones) 
of our sample and the Hamiltonian in this Gibbs ensemble is just given by 



E 

a=l,M 



H (a) = V H(a a ) . (38) 



Obviously for large M the correlation function can be defined as 

= SB=1 ^ ( *" K(t) (39) 

and in the limit M — > oo of this quantity is self-averaging, i.e. it will be the same in all the history of the 
system. It is evident that the M systems are independent so that the average of one Gibbsian systems is 
equivalent to the average of M usual systems and corresponds to repeat the same experiment (or computer 
run) M times. A similar procedure can be done for the response function. 
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5 Perturbing the system 

5.1 Global perturbation 

As usual we can rederive the probability distribution of the overlaps by perturbing the system. Let us 
consider for simplicity the effect of adding to the Hamiltonian H K (a) an extra perturbation given by 

A#« = h i< > ( 4 °) 

i=l,N,a=l,M 

where the variables h are Gaussian random variables with zero average and variance S. 
By integration by part we find that for given K 

IVF" = NM = ^ J dqPt{q){l q) ■ (41) 

Therefore we have the relation 

X (1)s £5^ 1 =^/* P ^ 1 -«) (42) 

for the susceptibility x 1 ■ The proof is identical to the standard one. On the other hand at e = stochastic 
stability implies that 

X « =(3 f dqP(q)(l - q) , (43) 

where P(q) is the the usual J average of the J dependent probability distribution. 

It is a trivial task to generalize the proof to the other susceptibilities. For example if we define 

AH^ = Yl hf k afa b k , (44) 

i=l,N,a=l,M,k=l,N,b=l,M 

where the variables h are Gaussian random variables with zero average and variance 5/(NM), we get 

x<2) = Ji^o =P I dqPt ^ 1 - q2) = l3 S ^ Pt((z)(1 - q2) ■ (45) 

In this way we can compute all the moments of the both the function Pt(q) and P(q) and they coincide. 
The two functions are equal. 

5.2 Local perturbation 

We can now repeat the same steps as before but locally. Let us consider for simplicity the effect of adding 
to the Hamiltonian Hk{<j) an extra perturbation given by 

o=l,M 

where the variables h are Gaussian random variables with zero average and variance 8. By integration by 
part we find that for given K 

IvF" = m = 6(3 J dqPl{q){l q) ■ (47) 



Therefore we have the relation 



for the susceptibility 
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Let us consider an aging system and add the perturbation AH^ at time t w . We have that 

Xi '(*) = ]™< Ai ^ )< = ]jj = ( 49 ) 

However in the limit small e in the dynamics does not depends on a so that we get 

x f\t) = Rdt w ,t) . (50) 

Assuming that 

Um X f ) W = xf ) , (51) 

we arrive to 

lim RiXt w ,t') = X f } = /? / dqPi{q){l-q) . (52) 
We can now copy mutatis mutandis the proof of the usual FDR. For example let us define 

AH^ = J2 h^ofoi, (53) 

a=l,M,b=l,M 

where the variables h are Gaussian random variables with zero average and variance 8/{M). The static 
susceptibility is given by 

a,&=l, iU 17 

If we assume for simplicity a Langevin type of evolution, the same steps of JTo|| give we have that 

X?\t) = 2j t t dt'C l (t',t)^A. (55) 

Assuming for simplicity the asymptotic scaling (this step is not necessary indeed the following formulae 
force its validity) we get 

lim X f\t) = 2 ( dCC Xi(C) , (56) 

t^oo J 

where we have defined: 

Mt)(C(t ,*)) = [-^r- ) (57) 

and 

Xi{C) = lim Xi(t)(C) . (58) 

t — *O0 

If also in this case the limits t — > oo and £ — ► may be exchanged we get 

lim x f\t) = X f ) = /? / dgP<(g)(l - <z 2 ) . (59) 
Generalizing the previous arguments we get 

x^d) = sj* dmi^ty- 1 )^^ = s J dec^x^c) = 

P J dqP t (q)(l - q s ) = sp J X^q*- 1 . (60) 

We thus arrive to the announced conclusion that 

Xi(C) = Xi{C) EE f dqP l (q)dq . (61) 
Jo 

This is the local relation among the static and the local fluctuation dissipation relations. 
A few remarks are in order 
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• If we take a sequence of larger and larger systems, the dynamical quantities converge to a well denned 
limit when the volume goes to infinity. Therefore also the local overlap distribution Pi(q) goes to a 
limit when the volume goes to infinity and depends only an the local environment (i.e. the ./ not 
too far from the point i). It is clear that all the problems due to the definition of the function Pj(q) 
infinite volume limit due to the difficulties of define the infinite volume limit of local observation fades 
away. It is quite possible that the construction we have presented in this paper may be useful in a 
rigorous approach. 

• The way in which we have coupled together our M clones is irrelevant. The only thing we need to 
introduce such a coupling is to correlate the states of one clones with the states of the the other clone. 



5.3 Bilocal Perturbations 

We are now ready for proving thermometric indistinguishability of the sites. We consider two far away sites 
i and k and we perturbation that depends both the spins at i and the spins at k. 
A typical example is 

A77 £ 2) = Yl h ai ' a *' a3 ' bl ' b2 a? 1 a% 2 (T? 3 (T b i 1 (T b i 2 . (62) 

a,! ,a 2 ,a 3 ,bi ,b 2 = l ,M 

where the variables h are Gaussian random variables with zero average and variance 5/M 4 . In the same 
way as before we get that 

(.«,.*) = p f dq , idqkP ( quqk){ l _ (g t ) S *(%) Sfe ) , (63) 



where P(qi, qk) is the probability distribution of qi and qk- If we compute the same quantity in the aging 
regime for very large time we get that the same quantity must be equal to the large time limit of 



dt 



r)C Si 3C Sk 



(64) 



Imposing that the two expressions can be equal for arbitrary and Sf. one recover the form of the P(qi, qk) 
in terms of Xi(t) and Xk(t). By imposing that P(qi, qk) is positive (i.e. it does not contain a derivative of 
a delta function) one finds that 

Xi(t) = X k (t) . (65) 

We finally obtain 



P(Qi,Qk) = 

Jo 



dx8(qi - q l (x))5{q k - qk{x)) , (66) 



where qi{x) is the inverse function of Xi(q). 

The first conditions is just thermometric indistinguishability of the sites during aging while the second 
condition has some interesting consequences that will be investigated elsewhere. 

We note that the probability distribution of the local overlap P%(q), being related to a dynamical 
quantity, must depend only on the local environment around the point i and therefore it must have a 
straightforward limit when the volume of the system goes to infinity (e.g it should be independent of the 
boundary conditions). It is remarkable that for far away points (i, k) the two probability distributions Pi{q) 
and Pk(q) are independent one from the other, but the joint probability distribution of qi and qk does not 
factorize has shown by eq. pq. 



6 Computing the local overlap distribution 

The formulae presented in the previous section are useful to define the local overlap and its distribution, 
but they are not handy as far as practical computations are concerned. The distribution Pi(q) will have 
rounded delta functions, for finite M. On the other hand the burden of the numerical computation increase 
violently with M. The best thermalization method (the parallel tempering) becomes more and more slow 
when the volume of the system (NM) increases. It is convenient to obtain this kind of information using 
different methods. 
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6.1 Computing the moments 



In many cases a the knowledge of the first few moments of the function Pi(q) is enough to compute this 
function. This approach is very useful in the case where we have some a priori information on the shape of 
this function. For example in the case of one step broken replica symmetry the knowledge of the first two 
moments and of m determines completely the function Pi(q). 

In other models, where the replica symmetry is broken in a continuos way it is possible that a good 
approximation is given by simple expressions like 

Pi(q) = 0(q - qf A )(l - 9 )- 1/2 g 4 (?) + m6(q - qf A ) (67) 

where Qi{q) is a low degree polynomial. However the validity of similar formulae may strongly change from 
model to model. 

The computation of the moments is rather simple. We consider the Gibbsian system with fixed M. If 
a and r are two equilibrium configurations of the model it is possible to prove that 

< n < T t) = i* M) = f d i i Mp M > (68) 

a=l,M J 

or more generally 



?<« . (69) 

where the indices a are arbitrary as soon as they are all different (obviously k < M). 

The proof of the previous relation can be obtained by computing the local susceptibilities like eq. (54). 
Let us consider the case k = 2. We could also define 

AtffW^Y 72 £ h^art, (70) 



The corresponding susceptibility should not depend on M and it is trivial equal to the previous defined 
susceptibility when M — > oo. In other words 

WRW) ■ (7i) 

does not depend on M, a and b as soon as a ^ b. 

This result can be checked in an explicit way in the two explicit cases discussed before eq. (2.4) . 

The computation of the moment of order k involves only the thermalization of a system with only kN 
sites and it can be done with not too much computational effort for not too large k. 



6.2 Introducing a state reservoir 

An other possibility, that is more interesting especially for analytic computations, consists in noting that 
the local overlap distribution depends only on the local environment. Therefore we can embed the local 
environment in a large system, and the local properties should not be related to that of a rest of the system 
as soon as the rest of the system has the same distribution of states of the original system. 

Let us consider a simple case, a spin glass model in three dimensions. Let us suppose that within the 
needed accuracy the knowledge of the couplings up to a distance R from the site i determines the function 
Pi(q). Let us consider a system of size L > 4R with the same probability distribution of the couplings, 
with the constraint however that there are two points i\ and ii such that the local environments of radius 
R around each of these two points (that are at distance greater than 2R) are the same and equal to that 
of the original system at the point i. 

In this case the same argument as before leads to the conclusions that if we take two different equilibrium 
configurations a and r, we have that 

(o-^cr^T^} = qf> , (72) 
where the bar denotes the average over the couplings that do not belong to the fixed environment. 
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A different possibility would be to take two copies of the local environment and to couple the spins at 
the boundary of each copy to other spins that stay on a Bethe lattice. If the state distribution of the Bethe 
lattice is the same as that of the three dimensional lattice, one should get identically results. 

In both case the main role of the system outside the local environment was to stay in different possible 
states, with the correct probability distribution. In the nutshell played the role of a state reservoir. As far 
as the detailed nature of this reservoir is not important we can consider the simplest possible model for it. 
A very convenient choice is the following: we model the interaction of the local environment with the rest 
of the world by introducing fields an interaction 



^ h k a k , (73) 



where the sum is done over the spins of the surface (S). The free energy corresponding to a given choice 
of the variables h is F[h] and the corresponding magnetization of the site i is m[h]. 

We should now modeling the ensemble of the h and of the w. We should introduce the distribution of 
the weight w a and for each state a we should give the set of the h. Let us consider the simplest case of a 
system with one step replica symmetry where we fix the value of m. In this case we can assume that the 



w are distributed according to the expression in eq. (18) . We can assume the variable h a have a given 



probability distribution, e.g. one of the simplest possibility is that 

ht = K + 5h% (74) 

where both h a and Sh^ are uncorrelated random variables with zero average and variance /i (0) and ftW. 
Different form of the probability distribution of the fields h can be used, however if the local environment 
is sufficient large, the result should not depend on the form of this probability distribution or from the 
variances and h&. 

One finally finds that the two parameters that identify the probability distribution, i.e. and q\ are 
given respectively 



E^m w " exp(-/3F[h a ])m[h a ] 



i = E a =i,A w « e M-PF[h a }Mh«}2 

qi E«=i,^ae X p(-/3F[^]) [ ] 

where A is a large number, the bar denotes the average over the magnetic field h and the weights w. Using 
the techniques of p2| in the one step replica breaking case the average over the w can be done and one 
obtains; 



T, a =i,A exp(-/3mF[fa"])m[fe° 
Ea=i,^exp(- / 3mF[/ l «]) 



l E a =i,A e M-P™F[h«])m[h«]2 
^ E a =i,A^M~PmF[h»}) [ ] 

Of course if our model is a Bethe lattice, the computations are quite simple and we reobtain the results 
of E2), i.e. a local version of the equations of pa]. 



7 Conclusions 

The main results of this note is the definition of a local probability distribution of the overlap, that depend 
on the site for fixed sample. The properties of this local probability distribution arc related to the local 
fluctuation dissipation relations, that automatical follows from the present formalism. The property of 
thcrmometric indistinguishability of the sites turns out to be a byproduct of our approach: during an aging 
regime all the sites are characterized by the same effective temperature during the aging regime. 
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The two times local correlations functions can be written 



CiXt w ,t)=Ci{x(t w ,t)) (77) 

where the function Ci(x) is the inverse of the function Xi(C) and can be obtained only by static measure- 
ments. The whole local dependence of off-equilibrium correlations and responses can be computed from 
static quantities. The only quantity that cannot be computed from equilibrium consideration is the global 
effective temperature, i.e. (3x(t w ,t). 

The physical picture that emerges is quite clear. The local overlap distribution can be defined by 
averaging over the ensemble of the states of the system and the introduction of the Gibbs ensemble is 
a technical tool to perform this average in a constructive way. During the dynamics the system locally 
explore different states of the system in a random way. It should not be surprising that the dynamical 
average is deeply connect to the static average on all possible different states. 

I am grateful to Jorge Kurchan, Marc Mezard and Federico Ricci-Tersenghi for illuminating discussions. 



Appendix 

It is interesting to consider what happens if we look to the model with Hamiltonian 



+ e 



= 1,M 



a,b=l,M;i,j=l,N 



isa,b ajb 



(78) 



in the limit M — > oo but with fixed N.. In this case the second part of the Hamiltonian coincide with that 
of the Sherrington-Kirkpatrick model. 

For small enough e one finds that the solution of the model is the replica symmetric one. It depends on 
a parameter q, that can be found by looking to the solution of the equations: 



q = ^ , (79) 

where the magnetization are obtained by considering the statistical average for one sample with Hamiltonian 

H(a) + hi<Ti (80) 

i 

and the magnetic fields are random independent identically distributed Gaussian variables with zero average 
and variance eq and the bar denotes the average over the fields h. 

However by increasing e this solution may becomes unstable (De Almeida Touless line) and at higher 
esp one has to look for replica broken solutions. A computation similar to the original one shows that 
replica broken phase happens when e 2 N = 0(1) and this result is at the origine of the condition e 2 N >> 1. 
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